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K  4  INTRODUCTION 
- — - 

^  Yin  recent  studies  of  phase  equilibria  in  refractory  compound  sy s terns j at 
ManLabs, I  explicit  methods  have  been  developed  for  predicting  the  vapor 
pressures  of  the  individual  component  elements  over  NaCl  type  metal  car¬ 


bides  , 


1,  2) 


These  methods  have  been  developed  and  applied  to  a  series  of 
refractory  binary  compounds  which  are  stable  over  a  wide  range  of  tempera¬ 
ture  and  composition.  Comparison  of  the  predicted  vapor  pressures  with 
experimental  measurements  of  congruent  vaporization  in  HfC  and  ZrC  yields 
excellent  results.  In  addition,  comparison  of  predicted  congruently  vapori¬ 
zing  composition  in  the  niobium-carbon  system  and  the  observed  congruency 
agree  quite  well. 

On  the  basis  of  this  experience,  it  seems  worth  while  to  extend  the  method, 
derived  for  binary  systems,  to  a  ternary  metal-metal-carbide  compound 
(Ta-Hf)C,  which  as  the  same  stoichiometry  and  crystal  structure  as  the 
binary  carbides  discussed  earlier. 

Under  these  circumstances,  a  theoretical  study  of  (Ta-Hf)C  was  under¬ 
taken  .  I  The  objectives  of  the  present  study  were  to 

a)  Apply  the  Schottky- Wagner  model'^  used  previously  to  describe 
binary  compounds  ,)^o  the  ternary  {Ta-Hf)C  in  order  to  generate 
equations  capable  of  predicting  the  vapor  pressures  of  tantalum, 
hafnium,  and  carbon  over  (Ta-Hf)C  as  a  function  of  temperature 


and  composition. 


..  /") 


f - 


b)  Relate  the  computed  vapor  pressure  curves  to  rates  of  vaporii^ation 
characteristic  of  specific  compositions  through  the  Langmiur 


c) 


and 


d) 


equation  ♦  [ 


Present  explicit  equations  for  the  vapor  pressure  and  rate  of 
vaporization  of  each  component  as  a  function  of  temperature  and 
composition. 

Carry  outlexplicit  numerical  evaluation  of  the  vapompressure  and 

rate  of  evaporation  of  each  component  at  2200^K,  2600^K,  2800  K 

A 

and  3000^K  and  at  carbon  levels  of  41%,  46%,  49%  and  50%  respectively 
for  comparison  with  experimentally  observed  result^^j  ^  ‘ 

Subsequent  portions  of  this  report  detail  the  methods  used  to  accomplish 
the  objectives  described  above. 
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II,  EXTENSION  OF  THE  SCHOTTKY- WAGNER  MODEL  TO 
tEKNARY  COMPOUND^  ■”""  " 

In  order  to  extend  the  treatment  of  Kaufman  and  Clougherty,^ 

performed  for  binary  compounds,  to  the  ternary  case,  we  consider 

a  system  A-B-C  containing  a  ternary  compound  having  the  O’ crystal 

structure.  The  composition  is  specified  by  setting  (  1-x-y)  equal  to 

the  atom  fraction  of  A,  x  equal  to  the  atom  fraction  of  B,  and  y  equal 

to  the  atom  fraction  of  C.  consider  a  case  where  two  sublattices 

exist,  the  A  and  B  atoms  occupying  one  sublattice,  while  the  C  atoms 

occupy  the  other  (i.e.  A  =  Ta,  B  =  Hf,  C  =  carbon).  If  stoichiometry 

corresponds  to  y  =  y^  and  the  total  number  of  sites,  filled  and  unfilled, 

is  N  then, 

N  .  =  Number  of  A  sites  =  (1-x-y)  (1-y,,)  (l-y)  N 
SA  O 

N  „  =  Number  of  B  sites  =  jc(l-y  )  (l-y)  ^N 

8  £5  O  => 

N  ^  =  Number  of  C  sites  = 
sC  'os 

while 

=  Number  of  A  atoms  =  (l-x-y)M 

N„  =  Number  of  B  atoms  =  xN 
B 

N^  =  Number  of  C  atoms  =  yN 
where  N  is  Avogardro^s  number. 

If  the  ratio  of  A  atoms  on  A  sites,  ^  atoms  on  B  sites, 

Nao»  is  equal  to  the  ratio  of  A  sites  to  B  sites  (ditto  for  the  B  atoms) 
then 

^AO  ""  atoms  on  B  sites  =  x(  1-x-y)  N  ( l-y)"  ^ 

N*  •  =  Number  of  A  atoms  on  A  sites  =  (I-x-y)^N(l-y) 

A  1 
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and 


Nbo  Number  of  B  atoms  on  A  sites  =  x(  l-x-y)N(  1-y)" 

2  •  1 

Nri  =  Number  of  B  atoms  on  B  sites  =  x  N{l-y) 

In  line  with  the  standard  state  convention  adopted  earlier,  the  free- 
energy  per  gram  atom  F^is  given  by  equation  (1) 

O'  ,  ,-r.  o  .  TT  o  .  V  O  4-  ^A+^A+  ^B+^B+/^C+^C+ 

F  =  (l-x-y)F^  +xFg  +yF^  +  (— )  AF  +_  + -^  +-^ 

.kT  inwp  (1) 

In  Eq*  1,  F^,  and  are  the  free  energies  of  pure  A,  B,  and  C 

at  the  temperature  in  question,  where 

^  0°K]  =  [  0^^]  =  0  at  one  atmosphere  (Z) 

is  the  reference  state.  Moreover,  AF^is  the  free  energy  of  formation  of 
the  ternary  compound  for  a  given  value  of  x  and  y^.  The  free  energies  of  A 
atoms  on  B  sites  and  B  atoms  on  A  sites  are  given  by  while 

Fa.j.i  free  energies  of  formation  for  A,  B,  and  C 

vacancies.  The  numbers  of  A,  B,  and  C  vacancies,  which  appear  in  Eq.  1 
are  given  as  follows: 

N^  ,  =  Number  of  vacant  A  sites  =  (l-x-y)(l-y  )N  (1-y)  ^l•(l-x-y)N 
j\  T  os 

^B+  ~  Number  of  vacant  B  sites  =  x(l-y^)N^(l-y)  ^ -xN 

N^ .  =  Number  of  vacant  C  sites  =  y  N  -  yN 
C+  ^o  s  ^ 

The  final  term  to  be  evaluated  in  Eq.  1  is  the  thermodynamic  probability 


factor  V/p  which  is  given  by, 

^sc' 


N  r,! 

sB 


^A+’  ^A1*^B0’  *  ^Bl’  ^Ao’ 
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Making  the  appropriate  substitutions  for  the  and  applying  Stirling's 
formula  and  substitution  into  Eq.  (1)  yields  for  the  case  y^  =  l/2 

(l-x-y)F^  ^ 

+(0. 5z  +  y-  1)  { l-y)“ \  ( 1-x-y)  F^^+  xF^^)  +  (0. 5z  -  y)Fj3^+x(  l-x-y)( l-y)"  V 

+RT(-z  In  0.5z  +  y  In  y  +  x  Inx  +  (1-x-y)  In  (1-x-y) 

+(0.5z  -  y)  In  (0,5z  -y)  +  {0.5z  +  y  -  1)  In  (0.5z  +  y  -  1)  )  (4) 

where  z  =  ratio  of  sites  to  atoms  =  N^/N,  W  =  F^q  +  ^BO’ 

of  formation  of  the  ternary  compound  AF  (x,  y^]  has  been  approximated  by  a 
linear  combination  of  the  free  energies  of  formation  of  stoichiometric  AC  and  BC 
which  represent  TaC  and  HfC.  Under  these  conditions  Eq.  (4)  reduces  to 
F^=  (1-x)  F^  +xFg°  +  2(Af/^.  +  x(AF^^  -  AF^^^_)  +  ^  F^.^)  +  V/  x(l-x) 

+RT  (x  In  X  +  ( 1-x)  In  (1-x)  )  (5) 

for  the  case  where  y  =  0  and  N  /N  =  2  which  corresponds  to  a  substitutional 

s 

solid  solution  of  A  in  B  with  the  C  sublattice  empty.  Eq.  (5)  is  recognizable 
as  being  the  regular  solution  approximation  for  substitutional  solid  solutions. 

In  the  other  limits  where  x  =  0  or  1-x-y  =  0  (i.e.  no  B  or  A  atoms  are 
present),  Eq.  (5)  reduces  to  the  equation  for  the  binary  compound  identical 
with  the  expression  derived  earlierl  ^ 

The  next  step  in  the  derivation  is  to  fix  the  composition  and  minimize  the 
free  energy  at  constant  temperature  and  pressure  by  letting  the  volume  (i.e. 
the  total  number  of  sites)  vary.  This  procedure  is  performed  by  setting  the 
derivative  of  F^(Eq*  4)  with  respect  to  z  at  constant  x,  y,  and  T  equal  to  0. 
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The  result  is 

(AFg"  -  Af/^))  +  (F^^+  =■( ‘-W '1^3+ -  +  F^^ 

=  -2RT  In  201  (6) 

where 

-2RT  In  2a'  =  RT  In  (0 , 5z)^  (0. 5z-y)”  ^  (0, 5z  +  y-  1)"  ^  (7) 

when  y  =  y  =  \/Z,  a  =  (z- l)/2z  =  0.5(N  -N)/N  .  Thus,  as  in  the  binary  case, 
o  s  s 

Cl'  is  the  fractional  number  of  vacant  sites  at  stoichiometry*  Substitution  of 

Eqs*  6  and  7  into  Eq,  (5),  i.e*  minimization  of  the  free  energy  with  respect 

to  volume,  at  constant  composition,  temperature,  and  pressure,  yields 

F°’=  (l-x-y)(F°-  F^^)  +x(F^-  F^^)  +  y(F^°  -  F^^)  +  x(  l-x-y)  ( 1-y)' 

+RT(x  Inx +  y  Iny  +  (l-x-y)  In  (1-x-y)  ln(0.5z-y) 

-(1-y)  In  (0.5z  +  y-1))  (8) 

Eqs.  6,7  and  8  completely  define  the  temperature  and  compositional 

dependence  of  the  free  energy  in  terms  of  the  parameters  F.  ,  F_  ,  F_,,  Oi 

Jdt  Ot 

and  W. 


For  the  case  of  carbide  type  O’ phase  (NaC  1  type  where  the  phase  boundary 

onthecarbon  rich  side  approaches  50%  carbon,  i.  e.  y^  =  l/2, )  we  can  approximate^ 
F^^  =  -RT  In  20!  =  -2AH^0°K]  -  RT  In  2  (9) 

where  ^  ^  i  n-  rr 

2AH  [0°K]  =  2AH^^[0°K]  +  x(  1-y)"  2AH^)  (10) 

^A+  =  -2Ah/^[0°K]  -  RT  In  2  -  ZAf/^,  [  T]  (11) 

and  - 

Fg^  =  -2AHg^[0°K]  -  RT  ln2  -  2AFg^[  T)  (12) 

Thus  in  cases  where  the  free  energies  of  formation  ofthe  compounds  AC  and 
BC  are  known,  at  stoichiometry,  asafuntion  of  temperature,  all  the  param¬ 
eters  required  are  specified  with  the  exception  of  W. 
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The  final  step  in  the  present  analysis  is  the  derivation  of  the  partial 
molar  free  energies.  For  the  case  of  a  binary  system  (i.e.  if  y  were  equal 


to  zero)  these  relations  are  well  known, 

8F 


and 


Fg=FMl-x)^ 


(13) 


(14) 


where  F*  and  F~  are  the  partial  molar  free  energies  of  A  and  B  and  x  is 
A  B 

the  atom  fraction  of  B  and  (l-x)  is  the  atom  fraction  of  A.  The  analogous 
expressions  for  the  ternary  case  are 


Fb  =  F  +  (l-x)  (||)y  -y^)^ 

.  xdf),,  +  (Uy)  (||)^ 


C  "  "  ■  ax  'y 
Since  Eqs. ,  6,  7,  and  8  yield 
O  ■*«-»  \  .  j  o 


(15) 

(16) 

(17) 


(|r)  =  .{F^  -  F^^)  +  (F^  -  Fg^.)  +  W{l-2x-y)  (l-y)*URT  Inx(l-x-y)"^ 


and 


(^) 

ay  X 


+0.5a(l.y)’^  <^B+*^A++2  -  AF^^.)  ) 

- (F^°-  Fa+)  +  W  x^(  1-y)"  0 . 5zx{  1-y )'  ^ F 

+Z(AF^^  -  ^F^q)  )  +RT  lny(0.5z  +  y-  1)  (1-x-y)"  ^(0.  5z-y)" ^ 


(18) 


B+"^A+ 


(19) 


Substitution  into  Eqs.  ,  15,  16,  and  17  yields  the  expression  for  the 
partial  molar  free  energies  as  follows: 
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+  W  (l-yf^  .0.5XZ  (l-yf^  +  2  (AF^^  -  AfJ^)) 


BC  — AC' 


+  RT  In  (1-x-y)  (0.5z  +  y- 1) 


^B"  ^B  =  *  W(l-x.y)^(l-y)*^  +  0.5(l-x-y)z(l-y)‘^(F^^-F 


B+  A  + 


+2(AFg^  -  +  RT  lnx(0.52  +  y-1)"^ 


F^^-  =  -  F^^^  +  RT  Iny  (0.  Sz-yf  ^  <22) 

Since  ^  ^ 

^A  -  Pa  I  yl / Pa°  <23) 

^B^-  t  y^  /Pb°  <24) 

e^nd  (T  (T 

-  F^=RT  lnpj3  [  X,  y]/p^  (25) 

where  [  x,  y]  ,  Pj^[  vl  »  yl  pressures  of  A,  B,  and 

C  respectively  over  the  alloys  and  p^,  P(^  are  the  corresponding 

vapor  pressures  of  the  pure  elements,  Eqs*  20  through  25  can  be  used  to 
compute  the  vapor  pressures  of  the  elements  over  the  alloys. 

In  particular  when  y  =  y^  =  0.5 

^A*^-  ^A°'-fA+*‘‘'*'=‘^  -2*<Fb+-Fa++2(AF®j.-AF;Jc)) 

+  RT  In  (0.5-x)  (l-2Q)a'^  (26) 


Fg°  =  -  F^^  +  4  V/  (0. 5  -  x)2  +  2  (0.  5  -  x)  (F^,  -  F^  ^  +2  (AF^^-AF^^)  ) 


B+  A  + 


BC  *AC 


+  RT  lnx(l  -  2  Ck)  a 


F^^-  F^  =  -  F^^^  +  RT  In  0.5  (1  -  2  a)  a"  ^ 


where  z  can  be  approximated  by  unity. 
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When  y  is  less  than  0.5, 
equal  to  4Q!^(l-y)^(l-2y)" 
duce  to 


z  =  2(l-y)  and  (0.5z  +  y  -  1)  is  approximately 
Under  these  conditions  Eqs.  19  through  21  re- 


*(l-y)-'(Fe^.F^^  +  2(iF^<.  -  AF®^)) 


+  RT  ln{l-x-y)(l-2y)  4' ^  a‘^( l-y)“^ 

=  -Fg^  +  W(  l.x-y)2(  l-y)“^  +( l-x-y)  ( l-y)"  ^Fg^ 
+  RT  lnx(l-2y)  4"^  ci^[Uy)~^ 


(29) 

(30) 


and 


F^  =  -  F^.^  +  RT  lny(l-2y)"^ 


(31) 


Z  2  •'1 

V/hen  y  is  greater  than  0.5,  z  s  2y,  and  (0.5z-y)  is  approximately  40!  y  (2y-l)" 


and 


-  F^°  =  -F^^+ Wx'^(l-y)”^  -xy(l-y)  ■'|Fb.-F^^  +2(AFg“^ 


+  RT  In(  l-x-y)  (2y- 1)  ^  (32) 

1^-Fg  =  -  Fg^  +  W(l-x-y)^(l-y)‘^  +  ( t-X-y)y(l-y)*^(F‘g+-F^^  +  2(AFg^^-AF^^)) 

+RT  lnx(2y-l)"^  (33) 

1^-  =  -  F^.^  +  RT  In  (2y-  1)  4“  ^  O'"  V  ^ 


As  pointed  out  previously,  these  equations  are  completely  explicit  only  if 
W  is  known.  Consequently  W  can  be  determined  from  limited  data  or  approxi¬ 
mated  by  W  =  0.  In  the  latter  case,  Eqs.  26  through  34  can  be  used  to  compute 
phase  equilibria,  vapor  pressures,  rates  of  free  evaporation,  or  compositions 
of  congruent  vaporization. 
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III.  APPLICATION  OF  MODEL  TO  Ta  -  Hf  -C  SYSTEM 

In  order  to  apply  the  equations  developed  in  Section  II  to  the  tantaluna- 

hafnium-carbon  system  in  explicit  terms,  it  is  necessary  to  assume  that  W  =  0. 

Moreover  since  the  NaCi  type  O'  phase  exists  for  values  of  y  (atomic  fraction 

of  carbon)  which  are  less  than  or  equal  to  0,5,  Eqs,  26  through  31  can  be  applied. 

In  this  system,  tantalum  is  element  A,  hafnium  is  element  B,  and  carbon  is 

element  C.  Table  I  contains  the  relevant  values  for  the  free  energy  of  formation 

of  stoichiometric  HfC  and  TaC.  These  values  are  the  same  as  the  values  em- 

(2) 

ployed  in  an  earlier  study  of  the  binary  carbides.  Consequently,  Eqs,.  (26) 
to  (28)  yield  the  following  expressions  for  y  =  0,5, 

RT£n  [x]  /  p®^^^  =  -2x(23,  100)  -  RTtna 

+RT£n{0.5-x)  (35) 

RTfn  =  "^H£+  +2(0. 5-x)  (23.  100)  -RTinQ! 

+RT£nx  (36) 

and 

RT£n  /P°C  ^  '^C+  "  ^TfnZ  a  (37) 

since  1  -  2a  is  nearly  equal  to  unity.  On  the  basis  of  Table  I  and  Eqs,  (9) 
and  (10) 

RTfnCtr:  -  35,900  -  x(l  -  y)"  ^  23,  100  cal/g.at,  (38) 
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Table  I 


Summary  of  Numerical  Values  Required  For 
Vapor  Pressure  Computations 


2AH^^^^[0°K]  z 

-35, 900  cal/mol 

0°K]  = 

-59,000  cal/mol 

zaf'^^cI  ’■1 

^Ta  + 

RT£n2 

^Hf+ 

cal/ mol 

cal/g.at. 

cal/g.at. 

cal/mol 

cal/g.at< 

2200 

-33,260 

+66, 126 

3034 

-54,200 

110, 165 

2600 

-32,780 

+65,095 

3585 

-53,  100 

108,515 

2800 

-32,540 

+64,579 

3861 

-52,  100 

107,240 

3000 

-32,300 

+64,063 

4137 

-51,200 

106,065 

source:  Ref(2)  pp{21)  and  (69) 

T°K 

1  O 

-log  Pc 

-log  P^a 

- 

n  O 

log  P 

(pressure  in  units  of  atmospheres) 

2200 

8.781 

11.206 

7.382 

2600 

6.  174 

8.346 

5,246 

2800 

5.  151 

7.219 

4,412 

3000 

4.266 

6.243 

3.692 

source:  C  • 

-  JANAF  Thermochemical  Tables 

Hf  and  Ta  -  S.A.R.  #2  AF33(657)8635  -  ’'Borides"  p  VIIf-6 
RTina=  -35,900  -  x{l-y)‘^  23,  100  (fromEqs.  (9),  (10)) 

T]  =  -  RT  jCn  2  Of  (from  Eq.  (9) ) 


^^Hf+  ■  ^Ta+  ”  ^^TaC^^  ""  23,  100  cal/g.at,  (fromEqs.  (11)  and  (12)) 
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V/hen  y  <  0.5,  Eqs.  29  to  31  yield 

RTiJnp^^^  [x.y  l/p^^a  =  "  ^Ta+  '  23,100  -  2RTi!na 

+  RT£n  (1-x-y)  (l-2y)4'^  (1-y)"^  (39) 

RT£np^^^[x.y)  =  +  (1-x-y)  (1-y)'^  23,  100  -  2RT£nC! 

+  RTin  X  (l-2y)4"Ml-y)‘^  (40) 

and 

RTUn  [x,  y]  / +  RT^n  y  (l-2y)"  ^  (41) 

Eqs.  35  through  41  define  the  vapor  pressures  of  tantalum,  hafnium, 
and  carbon  as  a  function  of  x,  y  and  T.  Numerical  calculation  of  these  quanti¬ 
ties  has  been  performed  for  four  carbon  levels,  y  =  0.50,  0.49,  0.46,  and  0.41 

at  temperatures  of  2200,2600,  2800  and  3000^K,  Graphical  repres  entations  of  these 

.  .  ■ 

results  are  shown  in  Figs,  1  to  16.  In  addition,  the  rates  of  Langmuir  evapora¬ 
tion  for  each  component  has  been  computed  as  follows: 

^  =  44.4  p^^^  [  x,y,T]  gms/cm^  sec  (42) 

V  44.4p^j^^[x,  y,  T]  T‘ ^ gms/cm^sec  (43) 

{>^  =  44. 4  p^  [  X,  y,  T]  T"  gms/cm^  sec  (44) 

where  ,  M„-,  and  are  the  atomic  masses  of  tantalum,  hafnium,  and 

1  a  rdt  ^ 

carbon  respectively. 


>!c  , 

These  temperatures  and  compositions  were  agreed  upon  at  a  meeting  between 
L,  Kaufman,  of  ManLabs  and  H.B.  Probst,  D.  Deadmore  and  I.  Zaplatynsky 
of  N. A . S . A . 
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log  vapor  pressure,  p,  (atmospheres)  ^  log  vaporization  rate  v  (gm. /cm.  sec.) 


Atomic  Percent  Hafnium 


Fi,.  1  ■  Vaporization  of  (Ta-Hf)  C  at  2200°K  and  50%  carbon. 
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log  vapor  pressure,  p,  (atmosphere)  ^  log  vaporization  rate  v  (gm./cm^  sec.) 


Fig-  2  -  Vaporization  of  (Ta-Hf)C  at  2200°K  and  49%  carbon. 
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log  vapor  pressure,  p,  (atmospheres)  ^  log  vaporization  rate  v  (gm./cm.  sec.) 


Fig,  3  -  Vaporization  of  (Ta-Hf)C  at  2200°K  and  46%  carbon. 
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log  vapor  pressure,  p,  (atmospheres)  ^  log  vaporization,  rate  v  {gm/cm^  sec) 


Atomic  Percent  Hafnium 


Fig.  4  -  Vaporization  of  (Ta-Hf)C  at  2200°K  and  41%  Carbon. 
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log  vapor  pressure,  p,  (atmospheres)  4  log  vaporization  rate  v  (gm/cm^sec.) 


Fig,  5  -  Vaporization  of  (Ta-Hf)C  at  2600  °K  and  50%  carbo 
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log  vapor  pressure,  p,  (atmospheres)  4  log  vaporization  rate  v  (gm/cm^  sec.) 


Atomic  Percent  Hafnium 

Fig. 6  -  Vaporization  of  (Ta-Hf)  C  at  2600°K  and 49%  carbon. 
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40 


50 


nium 

2600  °K  and  46%  carbon. 


log  vapor  pressure,  p,  (atmospheres)  ^  log  vaporization  rate  v  (gm/cm^sec.) 


L _ I _ I _ I _ I _ ^ _ LI 

0  10  2  0  30  40  50  60 


Atomic  Percent  Hafnium 

Fig.  8  -  Vaporization  of  (Ta-Hf)C  at  2600°K  and  41%  carbon. 
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g  vapor  pressure,  p,  (atmospher 


Fig.  10  -  Vaporization  of  (To-Hf)  C  at  2800°K  and  49% carbon. 
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T=2800“K 
46  ®/o  Carbon 


log  u 


log  vapor  pressure,  p,  (atmospheres)^  log  vaporization  rate,  o  (qm/crcr  sec) 


Fig,  12 -Vaporization  of  (Ta-Hf)  C  at  2800°K  and  41%  carbon. 
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T 


- r 

T=3000“K 
50  7«  Carbon 


/log  0 


log  vapor  pressure,  p,  (atmospheres) -r  log  vaporization  rate  i)  (gm/cm^  sec.) 


Fig.  15 .  Vaporization  of  (Ta-Hf)  C  at  3000°K  and  46%  carbon. 
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•ISIi 


T=3000‘’K 
41%  Carbon 


Vaporization  of  (Ta~Hf  )C  at  3000°K  and  41%  carbon 


The  total  rate  of  vaporization  ^  has  been  computed  by  summation  of  the  individual 


rates  i.e. 


^  =  ^Ta  ^  '^Hf 


Figs.  (1-16)  show  ^  [x,  y,  T]  in  addition  to  the  individual  pressures. 

Reference  to  Figs.  1-16  shows  that  for  most  of  the  composition  and  tem¬ 
perature  ranges  explored,  the  vapor  pressure  of  Ta  is  substantially  less  than 
Hf  or  C.  Thus,  all  of  the  compositions  will  lose  hafnium  and  carbon  relative  to 
tantalum.  Consequently  the  real  rate?  of  vaporization  should  change  with  time 
reflecting  changes  in  the  composition.  The  computed  values  of  0  correspond  to 
fixed  compositions. 

In  the  Hf-C  binary,  congruency  phenomena  are  predicted  and  have  been 
observed  (see  pp  21-38  of  reference  2).  These  conditions  correspond  to  unique 
cases  where  the  rates  of  carbon  and  hafnium  vaporization  are  such,  that  the 
ratio  of  carbon  to  hafnium  atoms  leaving  the  solid  is  the  same  as  the  ratio  of 
carbon  to  hafnium  atoms  in  the  soUd.  This  condition  is  represented  explicitly 

by  Eq.  (46  )• 


P&^c>yc»'^^  (46) 


’‘'in  Eq.  46,  the  vaporization  rates  are  in  moles/ cm  se 
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when  y  =  1-x,  i.e.  no  tantalum  is  present,  the  establishment  of  congruency 

yields  a  stable  composition.  This  situation  is  shown  graphically  in  Figs. 

(10)  and  (11)  of  ref  (2)  and  compared  with  experiment  in  Figs.  (18)  and  (19) 

of  ref  (2).  Thus  the  computed  rate  of  vaporization  for  a  congruency  situation 

should  be  stable  since  the  bulk  composition  does  not  change  with  time. 

Since  there  are  no  such  ternary  congruency  points  indicated  for 

Ta-Hf-C  alloys  (i.e.  no  solutions  for  Eq,  (46)  appear  except  when 

y  =  1-x  ),  it  is  expected  that  all  of  the  compositions  and  rates  of  vaporiza- 
c  c 

tion  will  be  time  dependent.  However,  it  is  possible  to  compute  *'quasi- 
congruency*’  conditions  corresponding  to  situations  where  the  ratio  of  Hf 
to  C  leaving  the  solid  is  in  the  same  ratio  as  hafnium  atoms  to  arbon  atoms 
in  the  solid.  This  computation  is  performed  by  applying  Eq.  (46)  to  Eqs. 

40  and  41  yielding. 

2RT£n  (l-2y^  )  =  RTj?n  ^h/^C 

+  +  3RT£nQ!-  (1-x  -y  )(l-y  )'^  23,100 

rll+  C  C  C 

+  RT£n  8(l-y^)^  (47) 

Eq.  (47)  can  be  employed  by  substituting  various  values  for  y^  ,  (i.e. 

y^  :=  •  ^c  "  0*495,  y^  =  0.49,  0.48  etc.)  and  obtaining  the  corresponding 

values  of  x  .  These  companion  values  of  x  and  y  can  then  be  substituted 
c  c  c 

into  Eqs.  39  to  45  to  compute  ^  f  x  ,  y  ,  T]  ,  Figs.  17  to  20  show  the 

c  c 

results  of  this  computation. 
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log  vaporization  rate,u  in  gm./cm^  sec. 


Atomic  Percent  Hafnium 


Computed  initial  rates  of  vaporization  at  2200 °K 

Fig.  17  -  Computed  Initial  Rates  of  Vaparization  at  2200®K. 


3 


T=2600“K 


in  gm./cm.  sec. 


T=3000“K 

@  Computed  on  basis  of  Hf-C  "conqruoncy " 
Vo  Carbon  indicated  iside  circle 
—  Computed  for  49%  C 


/  0.00%  To 

50.1 9%  Hf 
/\  49.81  %C 


/  Computed 


Atomic  Percent  Hafnium 


Computed  initial  rate  of  vaporization  at  3000° K 

Figp  20  -  Computed  Initial  Rates  of  Vaporization  at  3000^K, 


In  addition,  Figs.  17  to  20  repeat  the  i'  [x]  curves  corresponding  to 
49%  carbon  i.e.  y  =  0,49.  The  reason  for  including  the  49%  curves  is  as 
follows : 

Fig.  18  corresponding  to  2600^K  shows  the  results  of  vaporization 
data  obtained  by  Deadmore  at  NASA  on  a  series  of  stoichiometric  (50% carbon) 
Ta-Hf-C  alloys  at  2600°K.  These  data  are  the  results  of  2-4  hour  exposures 
representing  "average  vaporization"  rates.  The  Ta-C  results  at  2600*^K 
yield  5>  s  10*  gms/cm  sec,  which  corresponds  to  a  carbon  level  of  49%. 
Reference  to  Fig,  21  (which  is  the  companion  to  Figs,  14  and  15  of  ref,  2) 
shows  that  if  the  experimental  exposure  time  were  extended,  the  carbon  level 
should  decrease  resulting  in  a  further  reduction  in  i'.  This  reduction  of  i 
is  shown  in  Figs,  6-8.  However,  the  exposure  periods  in  the  TaC  case  was 
only  long  enough  to  produce  an  “average  ^  “  which  corresponds  to  an 
“instantaneous  ^  “  characteristic  of  a  49%C  alloy.  Consequently  as  a  first 
approximation  to  the  task  of  comparing  the  computed  values  of  ^  with  those 
observed,  49%C  is  taken  as  a  lower  limit  controlled  by  the  exposure  time. 
Thus  we  expect  that  the  “quasi-congruency**  v  corresponding  to  values  of 
Yc  equal  to  or  greater  than  0,49  can  be  reached  and  maintained,  but  that  the 
exposure  time  required  for  reaching  lower  values  of  y^  is  longer  than  the 
experimental  exposure  time. 

On  this  basis  the  *' computed**  ,  x,  y  curve  consists  of  two  branches: 

a)  C  =  i  [y^,  x^,  T]  for  y^^ 

h)  V  =  i  [x,  y=0.49,  T]  for  y  <  0.49 
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Vapor  pressure  of  Carbon 
over  pure  Graphite  = 


Calculated  compositional  dependence  of  the  vapor 
pressures  of  Ta  and  C  over  the  TaC^  phase  at 
2600 °K. 


CONCLUSIONS 


The  Schottky- Wagner  model  of  a  two -sublattice  compound  has  been  applied 
to  the  (Ta-Hf}C  ternary  compound  and  detailed  in  a  form  which  is  sixfficiently 
explicit  to  permit  computation  of  the  individual  vapor  pressures  of  the  elemental 
components  and  the  rates  of  Langmiur  evaporation  for  a  wide  range  of  compositions 
and  temperatures.  \  As  was  pointed  out  at  the  end  of  Section  III,jthe  computed 
rates  of  vaporization  represent  initial  or  instantaneous  rates  corresponding  to 
a  fixed  composition*  Th.e  exception  is  pure  HfC  which  exhibits  congruent 


vaporization*^ 

In  view  of  the  level  of  agreement  between  the  predictions  and  available 
data  (see  Fig*  18)  the  thermodynamic  description  of  this  ternary  compound 
detailed  in  Section  III  appears  quite  adequate,  fhe  comparison  of  predicted 
initial  rates  of  vaporization  and  observed  ratesJ(Fig.  18)^hows  fair  agreement, 
the  maximum  difference,  being  a  factor  of  5  at  2600  K.  In  addition,  the 
appearance  of  a  minimum  in  the  rate  of  evaporation  which  is  observed  near 
4TaC:  HfC  is  also  predicted  by  the  present  calculations.  Under  these  circum¬ 
stances  the  general  agreement  between  the  computed  results  and  the  experi¬ 


mental  behavior  appears  quite  satisfactory. 

The  present  computations  could  be  further  improved  by  converting  the 
equations  for  calculating  initial  rates  of  vaporizatio^,  Eqs.  (42>-(44),|to  inte¬ 
grated  rates  of  vaporization  which  might  be  compared  with  experimental  data 
on  a  more  reaUetic  basis.  [The  latter  problem  is  an  exercise  in  mathematics 
which  does  not  require  additional  thermodynamic  study. 
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For  example  consider  a  case  similar  to  the  binary  Ta-C  at  2600^K  (Fig.  25) 
where  the  pressure  of  metal  is  lower  than  the  pressure  of  carbon.  If  a  compo¬ 
sition  corresponding  to  y  =  0 . 5  (x  =  0)  is  prepared,  heated  to  2600^K  and  allowed 
to  vaporize,  then  the  initial  rates  of  vaporization  and  v  can  be  computed 
from  Eqs.  (35-44)  as  in  Fig.  (5).  However,  as  the  time  of  exposure  increases 
from  t  =  0  to  some  finite  time  t,  the  composition  will  change.  Since  congruency 
conditions  are  absent,  the  specimen  will  loose  carbon  and  y  will  decrease 
from  the  initial  value,  y^  =  0.5  at  t  =  0»  to  a  value  y[t]  which  is  less  than  0.5. 
As  this  compositional  change  occurs,  t  =  o]  will  decrease  (see  Figs.  6,  7,  8). 

In  order  to  compute  t]  and  v»[  t]  for  this  case,  the  changes  in  composition 
must  be  taken  into  account.  To  illustrate  the  indicated  procedure,  consider 
a  case  similar  to  the  TaC  situation  at  2600^K,  i.e.  x  =  o,  where  y  =  0.5  at 
t  =  0  •  Since  the  vapor  pressure  of  tantalum  is  much  lower  than  that  of  carbon 

.  If  and 

^Ta  number  of  grams  of  carbon  in  the  specimen  then 

y[  t]  =  t]  /  (w^°  ^Ta'  ^  ^  ^ 

where  w^^[  t]  =  constant  =  and  and  ^  are  the  atomic  weights  of 

carbon  and  tantalum  respectively;  when  t  =  0  ,  y[  t]  =0,5.  Eq,  (44)  defines 

^c[y]  =  -  ^.PJ.y2  =  44. 4  pjy]  =  y(  1.2y)^K^[  T]  (48) 

where 

K^[t]  =  exp[-F^^  /RT]  (49) 


initially,  let  us  assume  that 


L  tj 


W  cLS  cL 


ijLraL  d-pproximatiG 
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is  constant  at  a  given  temperature.  Substituting  (48)  into  (47)  yields 


^  =  -KjTl 


Integrating  yields 


°  lny/(l-y)  +  t 


hence 


(l-2y)/(l-y)  +  lny/(l-y)  =  -  t  K^/w^° 


Eqs.  (51)  and  (52)  define  y[  t]  and  the  carbon  weight  loss  (i.e.  the  total 
weight  loss  since  dWrp^/dT  is  assumed  equal  to  zero)  as  a  function  of  y  and 
time.  Combination  of  Eqs,  (51)  and  (52)  yields 


(Wc-Wc  ^  r  : 


(l-2y)/(l-y) 


Eqs.  51  to  54  can  be  simplified  by  considering  times  short  enough  so 
that  the  ln(l-y)/y  can  be  expanded  in  series  approximation.  This  procedure 


yields 


2  V-  1 


vjt]  r  2K^  (l-y)  (l-2y) 


tjt]  =  (2w^°K^t)^/2  (57) 

Thus  on  the  basis  of  the  present  simplified  expression,  the  rate  of  evapora¬ 
tion  is  found  to  vary  inversly  as  the  square  root  of  the  time,  A  more  complete 
solution  of  the  problem  would  require  relaxation  of  the  assumption  that 
dwTa/^^  =  o  • 


39 


REFERENCES 


1.  Kaufman,  E.  and  Clougherty,  EcV.,  "Thermodynamic  Factors 
Controlling  the  Stability  of  Solid  Phases  at  High  Temperatures  and 
Pressures".  Metallurgy  at  High  Temperature  and  High  Pressure, 
A.I.M.E.  Symposium,  Dallas  1963  to  be  published  by  Gordon  and 
Breech,  N.Y,  N.Y.  1964  Edited  by  K„  Gschneider,  N.  Parlee  and 
M.  Hepworth. 

2.  Kaufman,  L.  ,  Bernstein,  H.  and  Sarney,  A.,  "Thermodynamics  of 
Interstitial  Solid  Solutions  and  Refractory  Compounds".  A.S.D.-TR- 
61-445  pt  III  ManLabs  Report  November  1963  under  contract  No. 
AF33(657)-9826  with  Air  Force  Materials  Laboratory,  Dayton,  Ohio. 


40 


NASA-Langley,  1966  CR-458 


''The  aeronautical  and  space  activities  of  the  United  States  shall  he 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowU 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  ayid  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.” 

—National  Aeronautics  and  Space  Act  of  1958 


NASA  SCIENTIFIC  AND  TECHNICAL  PUBLICATIONS 


TECHNICAL  REPORTS:  ScientiEc  and  technical  information  considered 
important,  complete,  and  a  lasting  contribution  to  existing  knowledge. 

TECHNICAL  NOTES:  Information  less  broad  in  scope  but  nevertheless 
of  importance  as  a  contribution  to  existing  knowledge. 

TECHNICAL  MEMORAlNDUMS:  Information  receiving  limited  distri¬ 
bution  because  of  preliminary  data,  security  classification,  or  other  reasons. 

CONTRACTOR  REPORTS:  Technical  information  generated  in  con¬ 
nection  with  a  NASA  contraa  or  grant  and  released  under  NASA  auspices. 

TECHNICAL  TRANSLATIONS:  Information  published  in  a  foreign 
language  considered  to  merit  NASA  distribution  in  English. 

TECHNICAL  REPRINTS:  Information  derived  from  NASA  activities 
and  initially  published  in  the  form  of  journal  articles. 

SPECIAL  PUBLICATIONS:  Information  derived  from  or  of  value  to 
NASA  activities  but  not  necessarily  reporting  the  results  of  individual 
NASA-programmed  scientific  efforts.  Publications  include  conference 
proceedings,  monographs,  data  compilations,  handbooks,  sourcebooks, 
and  special  bibliographies. 


Details  on  the  availability  of  these  publications  may  be  obtained  from: 


SCIENTIFIC  AND  TECHNICAL  INFORMATION  DIVISION 

NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 

Washinston,  D.C  20546 


